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We study the thermal transport occurring in the system of solar captured dark matter (DM) 
and explore its impact on the DM indirect search signal. We particularly focus on the scenario of 
self-interacting DM (SIDM). The flows of energies in and out of the system are caused by solar 
captures via DM-nucleon and DM-DM scatterings, the energy dissipation via DM annihilation, and 
the heat exchange between DM and solar nuclei. We examine the DM temperature evolution and 
demonstrate that the DM temperature can be higher than the core temperature of the Sun if the 
DM-nucleon cross section is sufficiently small such that the energy flow due to DM self-interaction 
becomes relatively important. We argue that the correct DM temperature should be used for 
accurately predicting the DM annihilation rate, which is relevant to the DM indirect detection. 

PACS numbers: 


Dark matter (DM) composes about 25% of the energy 
density of in the universe and plays an important role 
in the structure formation. It was shown that if the 
galactic halo is constituted by weakly interacting mas¬ 
sive particles (WIMPs), there is a high possibility that 
these WIMPs are captured by the Sun [IHSIIIQ]. In gen¬ 
eral there are thermal energy flows between the captured 
DMs and the nuclei in the Sun. Microscopically, such 
flows are caused by particle scatterings. For collision¬ 
less cold dark matter (CCDM), the scatterings are only 
between the DM and solar nuclei. Regarding the huge 
difference in abundance between the two, it is reason¬ 
able to take the DM temperature to be identical to the 
core temperature of the Sun. However, the scenario of 
SIDM can change the picture dramatically. In such a 
scenario, the energy transports via DM self-capture and 
DM-nucleus scattering compete with each other. In par¬ 
ticular, if DM-nucleus interaction is much weaker than 
expected, the thermal exchange between DMs and nuclei 
would be much less efficient. As a result, the DM tem¬ 
perature can be distinct from the solar core temperature. 

It is worth mentioning that the DM abundance is not 
much affected with a suppressed DM-nucleon cross sec¬ 
tion a^p provided SIDM is considered m- The accu¬ 
mulated DM abundance is determined by the balance 
among the capture, the evaporation, and the annihila¬ 
tion rates. Earlier works on such processes only consider 
DM-nucleon interactions. The new effects from DM self¬ 
interaction are investigated recently The consid¬ 

eration of SIDM comes from the discrepancies between 
the numerical N-body simulations using the hypothe¬ 
sis of CCDM and the astrophysical observations on the 
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small structure of the universe [H]. The CCDM simula¬ 
tions m predict cuspy profiles in the center regions of 
galaxies, which conflict with flatten cores found in our 
Milky Way (MW) [2], other nearby dwarfs [T5], and 
low luminous galaxies |16[ I17j . There are additional puz¬ 
zles concerning the sizes of subhalos. The observed MW 
satellites are hosted by much less massive subhalos com¬ 
pared to sizes of the most massive subhalos arising from 
simulations. The absence of such massive DM subha¬ 
los is referred to as the “too big to fail” problem in the 
galaxy formation. The existence of SIDM (in particu¬ 
lar those involving the characteristic velocity-dependence 
cross sections [I8ll20j ) is one of the solutions to alleviate 
these inconsistencies. 

In this paper, we investigate the thermal energy trans¬ 
port between the trapped DM and the nuclei in the Sun, 
as well as the energy flows due to DM captures. With 
DM self-interactions taken into account, the DM number 
trapped in the Sun evolves according to 

^ = Cc + - CaNl , (1) 

while the energy flows in and out of the system of trapped 
DM in the Sun are governed by 


dt 


— Jc + {Jx + Js)Nx — JaN^ ■ 


( 2 ) 


The coefficients Cc,s are referred to as the DM capture 
rates due to DM-nucleon scattering cross section cr^p and 
DM-DM scattering cross section respectively. The 
coefficient Ca is related to DM annihilation rate in the 
Sun. In this work, we shall focus on the self-interaction 
dominant scenario with dCcCa [5]. The left hand 

side of Eq. ([^ is the total kinetic energy of the trapped 
DM with E^{t) the average kinetic energy of an indi¬ 
vidual DM. The factors Jc,s,a are thermal transport co¬ 
efficients corresponding to coefficients Cc, Cg, and Ca, 
respectively. The coefficient describes the heat ex¬ 
change between trapped DM and the nuclei in the Sun. 
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FIG. 1: The time evolutions of and for a^p — 10 cm^ (upper pannel) and a^p = 10 cm^ (lower panel) with 
^xx ~ 10“^® cm^ (left) and 10~^"^ cm^ (right). 


We do not concern the mass range where DM evapora¬ 
tion rate is significant. In this case, the DM abundance 
in the Sun is severely reduced so that the signal strength 
of DM annihilation is suppressed. 

Since we aim at studying the temperature evolution 
of the trapped DM, it is important to compare mean 
collision time between a pair of trapped DMs and that 
between a trapped DM and nucleus in the Sun. We make 
the comparison for the cases of spin-independent (SI) and 
spin-dependent (SD) DM-nuclei scattering cross sections, 
respectively. We first note that the mean collision time 
between two DMs in the Sun is 






(3) 


while the mean collision time between DM and nucleus 


in the Sun is 


Wo - 




tSI n 




(4) 


for the case of spin-independent cross section with Vq the 
solar volume, v the average velocity of trapped DMs, Ni 
the number of nucleus i in the Sun, and cr^A_ the spin- 
independent DM-nucleus scattering cross section. We ne¬ 
glect possible numerical factors on the right hand side 
of Eqs. ^ and (|^ for this order of magnitude estima¬ 
tion. The time scale for DMs in the Sun to reach 
thermal equilibrium can be estimated by the condition 
T®'^ ~ Wx('’'x'^)- solve this approximated equation, we 
first assume that is still far from the maximal 

value of and verify this assumption later. In this case, 
one can show that N^{t^^) = C'cT^ for Cf » ACcCa- 
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With this input, we obtain = ^JVqICc(t^^v. 

For further discussions, we shall take = 10 
GeV as a benchmark. For such a DM mass, we have 
~ iONna^p by assuming isospin invari¬ 
ant DM-nucleon couplings. Hence we may write r = 
= 401V//(T®p/crj^^A^^(T®'^). To simplify this re¬ 
lation further, we note that the average mass density 
of hydrogen in the Sun is roughly 1 g/cm^. In other 
words, Nh — 6 X 10®^ given the volume of the Sun 
approximately at 10^^ cm^. Using this value of Nh 
and the relation with Cc — 

5.4 X 10®®(cr®p/cm^)s“^ and v ~ 900 km/s, we obtain 

r ~ 10® We are interested in the parameter 

range that gives r < 1, i.e., < 10“^®. In this 

case, the trapped DMs reach to the thermal equilibrium 
among themselves before the energy exchange between 
the trapped DMs and the surrounding nuclei becoming 
efficient. Besides the condition r < 1, we need to ensure 
the ratio Cg/dCcCa = 1.9 x i® 

much greater than unity for the consistency of our ar¬ 
gument. The cross section combination = 

(10“'^®cm®, 10“®®cm®) is an example of satisfying both 
r < 1 and Cf/dCcCa ^ 1. With these parameters, we 
have = 4.5 x 10^® s, which is much shorter than 
the age of the Sun, r© « 10^^ s. For spin-dependent 
cross section, we found r ~ 1.8 x 10^ ^ and 

CH^CcCa = 10® X Since Cc — 

10®^ X (cr®p/cm®)s“®, the slightly different cross sec¬ 
tion combination = (5 x 10“^^cm®, 10“®®cm®) 

gives the same r®®* as the previous spin-independent case. 

The above discussions justify the thermal equilibrium 
state of DM in the early stage of capture for certain 
combinations of and We can write E^t) = 

skBTdt)/2 in Eq. ® where s is the degree of freedom of 
each DM. We recapitulate the meaning of various coeffi¬ 
cients in Eq. ([^. The coefficient Jc describes the energy 
flow due to DM capture caused by DM-nucleus scatter¬ 
ing, describes the energy exchange between DMs and 
nuclei, Jg is related to the energy flow due to DM capture 
caused by DM self-interaction, and Ja is related to the 
energy dissipation due to DM annihilation. All the quan¬ 
tities are positive except which depends on the differ¬ 
ence between the solar core temperature Tc and the DM 
temperature T^. In particular, < 0 when > Tc- 

When the Sun sweeps across the MW halo, there are 
collisions between DM and the solar nuclei. If the DM 
velocity is smaller than the solar escape velocity after the 
collision, the DM will be gravitationally bounded within 
the Sun. The velocity of an infalling DM at the shell 


^ In the subsequent discussions we shall only focus on spin- 
independent cross section since it is better constrained by direct 
detection experiments. We drop the superscript SI henceforth 
for simplicity. 


with radius r inside the Sun is tc = a/m® -I- ^^^^.(r). Here 
u is the DM velocity in the halo and n®gc(r) is the solar 
escape velocity at such shell. The kinetic energy of the 
trapped DM with an average over the energy loss in the 
capture process is given by 

- ( m^- wa V 2 I 2 

4 \ -I-TOA / 2 (m^-I- ttia)® ’ 

(5) 

where is the DM mass and toa is the mass of nuclei 
A. We have E = m;^c®s(,(r)/4 if it is due to DM-DM 
scattering, i.e., toa = rn^. We note that the collisions 
among DMs redistribute the DM kinetic energies such 
that the average DM kinetic energy becomes E^t) in 
Eq. ([^ when the thermal equilibrium is reached. The 
relation between E and E-^it) will be discussed later. The 
energy flow per shell volume is 


dJc 

w 



X 



(m^ — ttia)^ u® 
dm^ruA v%c{r). 


( 6 ) 


where ua and ct^a are the nuclei number density and the 
DM-nuclei cross section, respectively. The DM velocity 
in the halo is assumed to follow Maxwell-Boltzmann dis¬ 
tribution, i.e.. 


/(*) = , (7) 

V TT m^V XT] 

where a;® = 3(M/h)®/2 and ry® = 3{vq/v)‘^/2 with v « 
270 kms“^ the DM dispersion velocity in the halo and 
Vq = 220 kms“^ the relative velocity between the Sun 
and the MW. The density po is the DM local density® 
taken to be 0.3 GeV/cm®. The DM capture by Sun’s 
nuclei gives 


A 


(m® -I- m\) 
{m^ + Toa)® 



( 8 ) 


where 6 a is the number fraction of nuclei A in the Sun 
and {4>a) is the average gravitational potential square 
contributed by nuclei A. The quantity ^ is given by 


VesciRe) 3 /r. 

? = \/ o^QPO -=- VesciRe) - 


1.2 X 10®®GeVs“ 


Po 


0.3 GeV/cm® 


^270 km/s ^ 

(9) 


Similarly, the energy flow due to self-capture Jg can be 
derived by making the replacements toa -e and 


^ A slightly larger local DM density was derived in Ref. m- 
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riA —^ {rix number density in the Sun). 

We obtain 


7 ^ r^„ _ erf(?7) 7;esc(i?0)„,3 ^\2 

Js y 32 ^ 0 ^XX ^ - ^esci^&) {4’x) ■ (1*^) 

Since DMs are populated within a spherical region with 
a radius R ~ O.li?©, we can take (c^^) « (0x)^ where 
= 5.1 is the average gravitational potential of the 

DM. 

The energy of captured DM could be dissipated due to 
annihilation. The energy flow due to this process is 


^ J ATrr^nl(r)E^{t)dr 

(f 47rr^n^(r)dr)2 


( 11 ) 


with (av) the thermally-averaged DM annihilation cross 
section. Numerically, we have 


Ja(t) 


7.5x10"®^ GeVs"^ 


/ snix f ExU)\ 

(lOGev) ( GeV J 


( 12 ) 

where s is the degree of freedom of each DM. We shall 
take s = 3 for the numerical analysis later. 

Finally, the captured DMs continuously exchange en¬ 
ergies with the solar nuclei. The energy transport due to 
collisions between single DM and the surrounding nuclei 
in the Sun is given by [3] 


•^x 



ksjTQ - Ty) 
{nix + toa)^ 


X 



mp^ksT^ + 

m^niA 


1/2 


(13) 


where pc « 110 g/cm® is the core density of the Sun, 
/a is the mass fraction of nuclei A. We note that the 
above relation is derived by taking s = 3. We have nu¬ 
merically solved Eqs. Q and ([^. Both equations run 
from to = 10^® s, which is roughly the time scale for 
DMs in the Sun to thermalize according to our earlier 
arguments. As a first approximation, the initial condi¬ 
tion for E^{t) at to = 10^^ s is taken to be E in Eq. ([^ 
with and r averaged. We point out that E is the DM 
average kinetic energy before the thermalization. When 
DMs reach to the thermal equilibrium, they are popu¬ 
lated more closely to the solar core. Hence one expects 
Exito) > E hy the conservation of mechanical energy, 
since the gravitational potential energy of the thermal- 
ized DM is in general smaller than that of DM before 
thermalization. Nevertheless E^ito) and E differ within 
a factor of two, and in fact the DM temperature in the 
current epoch is independent of E^ito) as we shall argue 
later. 

The evolutions of DM number and DM temperature 
are presented in Eig. 0 . We take tr^p = 10 10 cm^ 

and a^x = 10“^®, 10“^^ cm^ for rrix = 5, 10, 50 GeV, 
respectively as the benchmark points. It is clearly seen 


that Nx at t ~ 10^® s is indeed far below the maximal 
This justihes our assumption earlier for deriving t®'^. We 
also see that the DM accumulation is enhanced by DM 
self interactions. In particular, the Nx enhancement is 
more significant for smaller nix and it is very sensitive to 
the ratio c^x/^’xp- 

Fig.g shows that DM temperature is higher than the 
temperature of solar core in the early stage for different 
DM masses. The DM temperature drops with time 
due to energy flow from DMs to solar nuclei. However, 
Tx is not always approaching to T^. Those cases with 
axp = 10“"^^ cm^ show such a behavior. Particularly, 
Tx/Tc 3.5 for Uxx = 10“^® cm^ and nix ~ GeV in 
the current epoch. We note that JoA^^ caused by the DM 
annihilation does not affect the DM temperature. Hence 
DM temperature evolution is controlled by parameters 
Jc, Cc, Js, Cs and Jx- However, as Nx accumulates, the 
effect by Jc and Cc become negligible. We thus have 

'^^^^Jx + Js-CsEx{t). (14) 


We observe that Ex {t) or (t) approaches to a constant 
value when the right hand side of Eq. (141 vanishes. This 
implies Jx —>■ 0, i.e., ^ Tc, if the self-interaction in¬ 

duced coefficients Js and Cs are negligible. On the other 
hand, does not approach to Tc if DM self-interaction 
cannot be neglected. We note that the gap between 
and Tc increases for heavier DM. This is because Jx scales 
as l/TO;^,for heavy DM. Thus, to balance Js — CsEx{t), 
Tx — Tc must be enhanced to compensate the 1 /m^ sup¬ 
pression in Jx- We argue that the DM temperature 
in the current epoch is hxed once the values for ct^xj 
Gxp and m-x are given. Essentially, Txit-o) is the solu¬ 
tion to Jx + Js — CsExito) = 0 with Jx the function of 
CTxp, nix and Exit), and both Js and Cs the functions of 
dxx- Glearly the initial condition Exito) does not affect 
Txito). Indeed we have varied Exito) within a reason¬ 
able range in our numerical analysis and found that Ex it) 
iTxit)) approaches to the same fixed point. 

We remark that the thermal equilibrium time scale 
which we choose as to, is obtained by a rough estimation. 
To study the effect from such an uncertainty, we also 
solve Eqs. Q and ([^ with to = 10^“^ s. We find that 
Txito) remains the same. 

It is important to note that Nx depends on the DM 
temperature Tx- This is due to the dependence of Nx 
on Ca- Hence the indirect detection signal could also 
be sensitive to Tx- We compare the annihilate rate 
Fa = CaNx/2 calculated with precise with that com¬ 
puted with the assumption Tx = Tc- The result is pre¬ 
sented in Fig. DM temperature affects both Nx and 
Ca in Fa. With Tx > Tc, we have CaiTx) < CaiTc) be¬ 
cause Ca oc Tx^^^ [2. As a result, we can easily see that 
NxiTx) > NxiTc) from Eq. ([T). As long as the enhance¬ 
ment on Nx can overwhelm tne suppression on Ca-, the 
annihilation rate Fa would be enhanced by adopting the 
precise DM temperature Tx- In the limit C^ ACcCa, 
Fa —>■ C‘ll2Ca when Nx in the Sun reaches to the maxi- 
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FIG. 2: Total annihilation rate Fa for a^j, = 10 cm^ and 10 cm^ with = 10 cm^ and 10 cm^. Dashed lines 
are those results calculated with = Tc. 


mum. One has rA('rx) > T^Tc) as a result. This is seen 
from the right panel of Fig. with a^p = 10“"^^ cm^ and 
= 10“^^ cm^. However, with = 10“^^ cm^ and 
(j^p kept the same, is still growing, i.e., ~ CcIq. 

Hence Fa — CaC'^t^l‘2., which implies rA(T^) < rA(T'c). 
We presented such a scenario in the right panel of Fig. 
as well. 

We have taken the Sun as an example to illustrate 
the possible temperature difference between DMs and 
their surrounding medium. We stress that the deriva¬ 
tion in this paper is generally applicable to other mas¬ 
sive celestial objects. The only required information is 
the DM dispersion velocity, the DM local density, and 
physical properties of the celestial object. In summary, 
we have derived and solved the thermal transport equa¬ 
tion for DMs trapped in the Sun for the first time. In the 


SIDM scenario, the DM temperature and the core tem¬ 
perature of the Sun could be different. We have shown 
that DM annihilation rate is sensitive to the DM temper¬ 
ature. Hence it is imperative to adopt the precise DM 
temperature to calculate the indirect DM event rates. 
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